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Abstract. The paucity of old isolated accreting neutron stars in ROSAT 
. observations is used to derive a lower limit on the mean velocity of neu- 

rS \ tron stars at birth. The secular evolution of the population is simulated 

j^ ■ following the paths of a statistical sample of stars for different values of 

the initial kick velocity, drawn from an isotropic Gaussian distribution 
with mean velocity < (V) < 550 kms~^. The spin-down, induced by 
dipole losses and the interaction with the ambient medium, is tracked 
together with the dynamical evolution in the Galactic potential, allowing 
for the determination of the fraction of stars which are, at present, in 
each of the four possible stages: Ejector, Propeller, Accretor, and Georo- 
tator. Taking from the ROSAT All Sky Survey an upper limit of ~ 10 
accreting neutron stars within ~ 140 pc from the Sun, we infer a lower 
bound for the mean kick velocity, (V) ^ 200 — 300 kms~^. The same 
conclusion is reached for both a constant {B ~ 10^^ G) and a magnetic 
field decaying exponentially with a timescale ~ 10^ yr. Present results, 
moreover, constrain the fraction of low velocity stars, which could have 
escaped pulsar statistics, to ^ 1%. 
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Isolated neutron stars (NSs) are expected to be as many as 10^-10^, a non- 
negligible fraction of the total stellar content of the Galaxy. The number of 
observed radio pulsars is now ~ 1,000. Since the pulsar lifetime is ~ 10^ yr, 
this implies that the bulk of the NS population, mainly formed of old objects, 
remains undetected as yet. Despite intensive searches at all wavelengths, only a 
few (putative) isolated NSs which are not radio pulsars (or soft 7 repeaters) have 
been recently discovered in the X-rays with ROSAT (Walter, Wolk & Neuhauser 
1996; Haberl et al. 1998; Neuhauser & Trumper 1999). The extreme X-ray to 
optical flux ratio (> 10^) makes the NS option rather robust, but the exact 
nature of their emission is still controversial. Up to now, two main possibilities 
have been suggested, either relatively young NSs radiating away their residual 
internal energy or much aged NSs accreting the interstellar medium (ISM) . Both 
options have advantages and drawbacks. Standard cooling atmosphere models 
fail to predict in a natural way the spectrum of the best studied object, RX 
J 1856-3754 (see Walter et al. this volume). Accretion models require instead a 
very low NS velocity relative to the ISM {v < 20 kms^^) in order to produce the 
luminosities inferred from ROSAT data (see Walter et al. this volume). We feel 
that a more thorough analysis of the statistical properties of NSs can be useful 
in providing indirect evidence in favor or against the accretion scenario. 

As discussed by Lipunov (1992), isolated NSs can be classified into four 
main types: Ejectors, Propellers, Accretors and Georotators. In Ejectors the 
relativistic outflowing momentum flux is always larger than the ram pressure 
of the surrounding material so they never accrete and are either active or dead 
pulsars, still spun down by dipole losses. In Propellers the incoming matter 
can penetrate down to the Alfven radius, Ra, but no further because of the 
centrifugal barrier, and stationary inflow can not occur, but the piling up of 
the material at the Alfven radius may give rise to (supposedly short) episodes 
of accretion (Treves, Colpi &: Lipunov 1993; Popov 1994). Steady accretion is 
also impossible in Georotators where (similarly to the Earth) the Alfven radius 
exceeds the accretion radius, so that magnetic pressure dominates everywhere 
over the gravitational pull. It is the combination of the star period, magnetic 
field and velocity that decides which type a given isolated NS belongs to and, 
since both P, B and V change during the star evolution, a NS can go through 
different stages in its lifetime. 

While the dynamical evolution of NSs in the Galactic potential was studied 
by several authors (sse e.g. Madau & Blaes 1994; Zane et al. 1995), little 
attention was paid to the NSs magneto-rotational evolution. Recently, this 
issue was discussed in some detail by Livio, Xu & Frank (1998) and Colpi et al. 
(1998). Goal of this investigation is to consider these two issues simultaneously, 
coupling the dynamical and the magneto-rotational evolution for the isolated 
NS population. 

The possibility that the low-velocity tail is underpopulated with respect to 
what was previously assumed should be seriously taken into account. It is our 
aim to revise the estimates on the number of old accreting neutron stars in the 
Galaxy in the light of these new data, in the attempt to reconcile theoretical 
predictions with present ROSAT limits (Neiihauser & Triimper 1999). 
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2. The Model 

In this section we summarize the main hypothesis introduced to track the evo- 
lution of single stars and describe shortly the technique used to explore their 
statistical properties, referring to Popov & Prokhorov (1998) for details on spa- 
tial evolution calculations and to Konenkov & Popov (1997) and Lipunov & 
Popov (1995) for details of magneto-rotational evolution. 

2.1. Dynamical evolution 

The dynamical evolution of each single star in the Galactic potential (taken in 
the form proposed by Miyamoto &: Nagai 1975) is followed solving its equations 
of motion. 

The period evolution depends on both the star velocity and the local density 
of the interstellar medium, any attempt to investigate the statistical properties 
of the NS population should incorporate a detailed model of the ISM geogra- 
phy. Unfortunately the distribution of molecular and atomic hydrogen in the 
Galaxy is highly inhomogeneous. Here we use the analytical distributions from 
Bochkarev (1992) and Zane et al. (1995) for the hydrogen density n{R,Z). 
Within a region of ~ 140 pc around the Sun, the ISM is underdense, and we 
take n = 0.07 cm^^. 

In our model we assume that the NS birthrate is constant in time and 
proportional in magnitude to the square of the local gas density. 

Neutron stars at birth have a circular velocity determined by the Galactic 
potential. Superposed to this ordered motion a kick velocity is imparted in a 
random direction. We use here an isotropic Gaussian distribution (relative to 
the local circular speed) with dispersion ay, simply as a mean to model the 
true pulsar distribution at birth (see e.g. Cordes & Chernoff 1998). The mean 
velocity {V) = {8/7:)^''^av is varied in the interval 0-550 km s^^. 

2.2. Accretion physics and period evolution 

The accretion rate was calculated according to the Bondi formula 

where rup is proton mass, the sound speed Vs is always 10 kms~^ and fio = 
(F^ -|- Vg)^''^ in units of 10 kms~^. M and R denote the NS mass and radius, 
which we take equal to 1.4 Mq and 10 km, respectively, for all stars. 

All neutron stars are assumed to be born with a period -P(O) = 0.02 s, and 
a magnetic moment either ^30 = 1 or n^o = 0.5, where ^30 = fi/W'^^Gcnr^. 

In the ejector phase the energy losses are due to magnetic dipole radiation. 
When the gravitational energy density of the incoming interstellar gas exceeds 
the outward momentum flux at the accretion radius, Rac — 2GM/v'^, matter 
starts to fall in. This happens when the period reaches the critical value 

PEiE^P)^ 10 ^li'n-'/^vli's. (2) 

When P > Pe{E -^ P) the NS is in the propeller phase, rotational energy is 
lost and the period keeps increasing at a rate taken from Shakura (1975). 
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Figure 1. Fractions of NSs in the different stages vs. the mean kick 
velocity for /i3o = 0.5 (open circles) and /X30 = 1 (filled circles); typical 
statistical uncertainty for ejectors and accretors is ~ 1-2%. 



As the star moves through the inhomogeneous ISM a transition from the 
propeller back to the ejector phase may occur if the period attains the critical 
value 



PE{P^E)^'ilx^^^'v 



f/^,6/^.-2/7 



10 



n 



s . 



(3) 

Note that the transitions P ^ E and E ^f P are not symmetric as first discussed 
by Shvartsman in the early '70s. 

Accretion onto the star surface occurs when the corotation radius Rco = 
{GM P'^ /A'k'^Y'^ becomes larger than the Alfven radius (and Ra < Rac: see 
below). This implies that braking torques have increased the period up to 



Pa{P ^ A) ~420/i3(,'n 



6/7„-3/7,9/7 



-'10 



S . 



(4) 



As soon as the NS enters the accretor phase, torques produced by stochastic 
angular momentum exchanges in the ISM slow down the star rotation at the 
equilibrium period 



P. 



eq 



2.6 X 10-' V 



3 ,,-2/3, 2/3^-2/3,13/3 



(t)10 /^30 ^ ^10 



(5) 
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Figure 2. Number of accreting NSs in the Solar vicinity for a constant 
(/"30 = 1) filled circles; /iso = 0.5, opened circles) and decaying field 
{td = 2.2 X 10^ yr, triangles). 



where vu\ the turbulent velocity of the ISM (Lipunov & Popov 1995; Konenkov 
& Popov 1997). 

At the very low accretion rates expected for fast, isolated NSs, it could 
be that the Alfven radius is larger than the accretion radius. The condition 
Ra < Rac translates into a limit for the star velocity 



V <410n^/^V3o^/^ kms"^ 



(6) 



3. Results and discussion 

3.1. The NS census for a non decaying field 

We consider two representative values for the (costant) magnetic dipole moment, 
^30 = 0.5 and ^30 = 1. The present fraction of NSs in the Ejector and Accretor 
stages as a function of the mean kick velocity is shown in figure 1. 
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Figure 3. Fractions of NSs in the different stages vs. the average 
kick velocity for a decaying field with an e-folding time td = 2.2 x 10^ 
yrs (open circles) and t^ = 1.1 x 10^ yrs (filled circles). 



Here, and in the following the total number of Galactic NSs was assumed to 
be 10^. A total number ~ 10^ appears to be consistent with the nucleosynthesis 
and chemical evolution of the Galaxy, while 10^ is derived from radio pulsars 
observations. It is uncertain if all NSs experience an active radio pulsar phase, 
due to low initial magnetic fields or long periods, or to the fall-back in the 
aftermath of the supernova explosion. There is a serious possibility that the 
total number of NSs derived from radio pulsar statistics is only a lower limit. 

In order to compare the expected number of accreting ONSs with the 
ROSAT All Sky Survey (RASS) results, we evaluated the number of those ONSs, 
within 140 pc from the Sun, producing an unabsorbed flux of 10~^^ erg cm~^ s~^ 
or higher at energies ~ 100 eV. The results are illustrated in figure 2. The main 
point is that for mean velocities below 200 km s~^ the number of ONSs with a 
flux above the RASS detection limit would exceed 10. Most recent analysis on 
the number of isolated NSs in the RASS (Neiihauser & Triimper 1999) indicate 
that the upper limit is below 10. 

An important aspect is that our results exclude the possible presence of a 
consistent low-velocity population at birth, which exceeds that contained in the 
gaussian with {V) > 200 km s~^. 
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Figure 4. Ejector time vs. bottom magnetic momentum for two 
values of initial magnetic momentum and different decay timescales, 



td = 10^, 10^, and 10 



10 



yr. 



3.2. The NS census for a decaying field 

We refer here only to a very simplified picture of the field decay in which B{t) = 
B{0)exp{—t/t(i). Calculations have been performed for td = 1.1 x 10^ yr, t^ = 
2.2 X 10^ yr and /U3o(0) = 1. Results are shown in figure 3. 

For some values of td and bottom field most of NS can stay at the Ejector 
stage, and the number of accretors and Propellers would not be increased. We 
show this analytical estimates graphically in figure 4, where the Ejector time, 
Te, is plotted vs. bottom magnetic momentum for constant velocity and ISM 
density (n = 1 cm~^, v = 10 km/s), different td and two values of the initial 
magnetic momentum, 10^'^ and 10^^ G cm^ (see Popov &: Prokhorov 1999). 

Summarizing, we can conclude that, although both the initial distribution 
and the subsequent evolution of the magnetic field strongly influences the NS 
census and should be accounted for, the lower bound on the average kick derived 
from ROSAT surveys is not very sensitive to B, at least for not too extreme 
values of td and /u(0), within this model. 
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4. Conclusions 

In this paper we have investigated how the present distribution of neutron stars 
in the different stages (Ejector, Propeller, Accretor and Georotator) depends on 
the star mean velocity at birth. On the basis of a total of ~ 10^ NSs, the fraction 
of Accretors was used to estimate the number of sources within 140 pc from the 
Sun which should have been detected by ROSAT. Most recent analysis of ROSAT 
data indicate that no more than ~ 10 non-optically identified sources can be 
accreting ONSs. This implies that the average velocity of the NS population at 
birth has to exceed ~ 200 kms^^, a figure which is consistent with those derived 
from radio pulsars statistics. We have found that this lower limit on the mean 
kick velocity is substantially the same either for a constant or a decaying i?-field, 
unless the decay timescale is shorter than ~ 10^ yr. Since observable accretion- 
powered ONSs are slow objects, our results exclude also the possibility that the 
present velocity distribution of NSs is richer in low-velocity objects with respect 
to a Maxwellian. The paucity of accreting ONSs seem therefore to lend further 
support in favor of neutron stars as very fast objects. 
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